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AN EXPERIMENTAL INVESTIGATION 
OF WALL BOUNDARY-LAYER TRANSITION REYNOLDS NUMBERS 

IN AN EXPANSION TUBE 

By K. James Weilmuenster 
Langley Research Center 

SUMMARY 

Experimental measurements of boundary -layer transition in an expansion -tube 
test-gas flow are presented along with radial distributions of pitot pressure. An integral 
method for calculating constant Reynolds number lines for an expansion-tube flow is 
introduced. Comparison of experimental data and constant Reynolds number calculations 
has shown that for given conditions, wall boundary -layer transition occurs at a constant 
Reynolds number in an expansion -tube flow. Operating conditions in the expansion tube 
were chosen so that the effects of test-gas nonequilibrium on boundary-layer transition 
could be studied. 


INTRODUCTION 

The development of the expansion tube as a hypervelocity facility has entailed the 
study of several effects which cause the resulting flow parameters to differ, from their 
ideal values. One such phenomenon, reported in reference 1, is that under certain 
operating conditions the test-gas flow experiences a rather sudden drop or dip in pitot 
pressure. It was at first thought that this dip in pitot pressure could be attributed to the 
secondary diaphragm since the energy loss from the flow was of the same order as that 
required to rupture the secondary diaphragm. However, parametric studies involving 
the strength and mass of the secondary diaphragm show no effect on the dip. Mixing of 
the acceleration and test gases offered another possible explanation for the pitot -pressure 
drop. 

Mir els' (ref. 2) theory of shock attenuation shows that in order to maintain a mass 
balance, a great deal of the shocked acceleration gas is fed into the boundary layer which 
lies between the tube wall and the test gas. Any mixing of the shocked acceleration gas 
and the test gas would produce pitot pressures which are lower than those predicted for 
the test gas. Thus, the reduction in the pitot pressure might indicate the diffusion of 
acceleration gas contained in the wall boundary layer into the free stream. Extensive 



tests to determine whether the test gas was indeed contaminated by the shocked acceler- 
ation gas were negative. 

As suggested by the authors of reference 1, the transition from a laminar boundary 
layer to a turbulent boundary layer on the tube wall in the test gas along with the more 
rapid boundary -layer growth and subsequent closure is proposed as an explanation for 
the drop in pitot pressure. Thus, the ability to track and predict transition would be 
important as a parameter in determining expansion-tube test time and operating condi- 
tions. This paper is concerned with the experimental tracking of wall boundary -layer 
transition through the expansion tube and the calculation of transition Reynolds numbers. 
Data for this study were taken in the 9. 52 5 -cm -diameter Langley pilot model expansion 
tube which is described in reference 1. 


SYMBOLS 

Primary measurements have been made in U.S. Customary Units and converted to 
SI Units. 

a speed of sound, meters/second 

b integrated Reynolds number parameter defined by equation (2) 

p pressure, newtons/meter^ 

P t pitot pressure, newtons/meter^ 

R Reynolds number 

Rj general integrated Reynolds number through the expansion fan 

Rf 0 integrated Reynolds number of the reference particle due to traversing the 

expansion fan 

R^j unit Reynolds number in region i 

R^ total integrated Reynolds number 

R (2),fan = R © + % 
r tube radius 
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Si 


shock propagating into region i 


t time, seconds 

U velocity, meters/second 

x distance, meters 

y distance measured from tube wall 

Subscripts: 

<t center line 

ent expansion -fan entrance conditions 

ex expansion -fan exit conditions 

f regions of flow affected by an unsteady expansion 

i regions in flow at different thermodynamic states denoted by (D . . . 

o reference conditions 

p particle property 

tr reference to transition 

Primes denote coordinates referenced to apex of centered expansion fan. 

EXPANSION-TUBE OPERATION 

In the operating cycle of an expansion tube, a very short duration steady flow is 
produced at hypersonic conditions. An idealized schematic of the expansion -tube flow 
cycle is shown in figure 1. The principles of operation are discussed in detail in 
reference 3. 

Briefly, the facility consists of a high-pressure driver section which is charged 
with either hydrogen or helium, an intermediate section which is charged with a test 
gas of interest, and an acceleration section which is normally charged with helium at a 
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low density. A heavy steel diaphragm separates the driver and intermediate sections 
and a thin plastic diaphragm separates the intermediate and acceleration sections. A 
shock-tube flow is initiated when the primary diaphragm is ruptured. The incident 
shock then ruptures the secondary diaphragm and induces another shock-tube -like flow, 
where the driven gas of the first flow becomes the driver gas for the second flow. Gas 
processed by the shock in the acceleration section is separated by an interface from the 
test gas (region (5)). A history of the processed test gas is indicated by a particle path 
originating in the intermediate section. This path represents a picture of the ideal, 
inviscid expansion-tube cycle. The inclusion of viscous and real-gas effects can signifi- 
cantly alter the previous description primarily through shock attenuation as described by 
Mirels (refs. 2 and 4) and nonequilibrium effects reported by Haggard (ref. 5). 

The experimental data that are presented in this paper were taken from a series of 
runs in the Langley pilot model expansion tube. The driver, test, and acceleration gases 
as well as the initial pressure in the driver and acceleration sections of the expansion 
tube remained the same throughout the tests. The driver gas was H 2 at a pressure of 
8.9 MN/m^ and the acceleration gas was He at a pressure of 0.27 kN/m^. The only 
parameter varied was the initial pressure of the test gas which was air. This change in 
pressure had the effect of causing a significant change in the unit Reynolds number in 
regions ( 0 ) and © although only varied from 3698 m/sec to 4622 m/sec. 

Flow parameters measured at the test section were made at different axial loca- 
tions downstream of the second diaphragm by moving the diaphragm station relative to 
the fixed test section. 


PITOT-PRESSURE DIP PHENOMENON 

Dips in the pitot pressure in region (D were first noted in reference 1 on center - 
line pitot-pressure records taken in the test section as illustrated by the lower oscillo- 
gram trace in figure 2. 

The small initial pressure rise is caused by the arrival of the shock S (g), 
followed by the region (Q$) (acceleration gas), while the second rise is caused by the 
arrival of the interface, with the subsequent region (D (test gas) following. The pitot - 
pressure level in region © should persist for approximately 200 jusec at a constant 
level and then rise as the unsteady expansion arrives at the test section. As seen in 
figure 2, the region d) pitot pressure remains constant for approximately 40 psec, drops 
sharply, and then begins to rise again. The upper trace in figure 2 is the wall static 
pressure taken at a station 26.5 cm upstream of the test section. When plotted on an 
x,t diagram, the high-frequency pressure fluctuations seen on the wall static pressure 
correlate well with the initiation of the dip in pitot pressure measured on the center line. 
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The high-frequency pressure fluctuations in the wall static pressure are indicative of the 
situation found in a turbulent boundary layer. However, as shown in figure 3, expansion- 
tube operating conditions have been found so that there is no apparent dip in the center - 
line pitot pressure whereas the wall static pressure still indicated the onset of the high- 
frequency pressure changes. Thin -film heat -transfer gages were placed in the tube 
wall to be used as further indication of transition. There was a good correlation between 
the time of onset of wall pressure fluctuations in figure 3 and the time of boundary -layer 
transition as indicated on wall -temperature records in figure 4. The time at which 
transition occurs in figures 3 and 4 is measured relative to the time of arrival of the 
incident shock as noted in the respective figures. Further investigation showed, as 
indicated in figure 5, that at the same conditions as the data presented in figures 3 and 4, 
pitot-pressure samples taken at one-half the tube radius above and below the center line 
indicated a presence of the pitot-pressure dip, an indication that the disturbance causing 
the dip moved from the tube walls toward the center line. Also, an order -of -magnitude 
reduction in the region unit Reynolds number used for the data shown in figures 3 to 5 
was found to eliminate pitot-pressure dip as shown in figure 6. The exact nature of the 
disturbance seen on all the records in figure 6, approximately 800 psec after the arrival 
of the incident shock, is not clear. It could be related to boundary -layer transition or to 
upstream moving waves associated with the normal operating cycle of the expansion tube. 
Similar data were reported by Jones and Moore. (See ref. 1.) 

Since the pitot-pressure dip moves from the tube walls to the center line and 
appears to be associated with indications of boundary -layer transition, it was believed 
that the pitot -pressure dip is triggered by the arrival of the outer edge of a rapidly 
growing turbulent boundary layer at the radial location of a pitot -pressure probe because 
of the reduction of total pressure and velocity in the boundary layer. The same phenom- 
enon would occur if the boundary layer were laminar since, for a given tube diameter, 
boundary -layer closure would occur after a sufficiently long time. However, this con- 
dition is unlikely to occur because of the short times involved and the character of the 
boundary layer in an expansion-tube -like flow as described by Trimpi and Cohen. (See 
ref. 6.) 

In an attempt to understand the dip phenomenon better, a five -probe pitot -pressure 
rake was designed and tested in the Langley pilot model expansion tube. The rake was 
designed so that the pressure transducers were located near the leading edge of the rake; 
thus, the response time of the probe was minimized. This rake sampled the pitot pres- 
sure on the center line and in four equal increments to a position y/r = 0.275. Thus, 
a radial sampling of the expansion -tube flow was obtained. Some of the results from 
these tests are shown in figures 7 to 9 for three different axial locations. Because the 
flow can be sampled at only one axial location per run, the data for each axial location 
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correspond to a different run. Efforts were made to maintain the same flow properties 
for each test run. Each of the pitot-pressure profiles is normalized with respect to the 
center-line pitot pressure. The oscillations in the pitot pressure have been eliminated, 
and the resulting smoothed curves have been plotted at the top of each figure. 

Several difficulties were encountered in the acquisition of the pitot-pressure data, 
the most significant being the acceleration sensitivity of the pressure transducers used 
in the rake. This acceleration sensitivity caused an estimated 20- to 30 -percent uncer- 
tainty in the data. Therefore, the data as shown in figures 7 to 9 are only intended to 
show possible trends in the development of the expansion -tube test-gas flow. It may be 
noted that the trend is for the dip to occur earlier in time relative to the interface as the 
point of data sampling is moved downstream of the secondary diaphragm, and for the 
wall boundary layer to grow at a more rapid rate in the vicinity of the dip. Also, the 
pitot pressure indicates that the tube is almost filled by the boundary layer when the dip 
occurs; this condition might explain why the effects of transition are seen on the center 
line such a short time after transition is detected at the wall. 

The following sections of the paper are based on the hypothesis that the wall static- 
pressure and heat-transfer records provide adequate definition of wall boundary -layer 
transition. Next, an effort will be made to determine the proper transition Reynolds 
number to use and to determine whether transition occurs at a constant Reynolds number 
or is a function of unit Reynolds number. 

RESULTS AND DISCUSSION 

In a flow as complex as that found in the expansion tube, an, a priori analysis to 
determine transition Reynolds numbers would be, at best, exceedingly difficult and no 
attempt will be made here to carry over steady -flow criteria to unsteady flow. As an 
alternative, a plausible method of calculating the pertinent Reynolds number for the 
experimental data is presented. 

When working with boundary -layer transition in a shock tube, Hartunian et al. 

(ref. 7) based the transition Reynolds number on distance a particle travels in the shocked 
gas before it reaches the transition point, which is equivalent to the path of the transition 
point being parallel to the path of the shock. Although there was some scatter in the 
experimental data, the use of this Reynolds number correlated his data very well. One 
might then expect the locus of points representing transition in the expansion -tube test 
gas to be at some fixed or nearly fixed distance from the contact surface. A plot of 
experimental transition points on an x,t diagram is shown in figure 10. As can be seen, 
the data hardly satisfy the previous assumption. The path of transition moves toward 
the contact surface or interface and eventually intersects it. 
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If a return is made to the idea of basing the Reynolds number on the distance a 
particle has traveled since being set in motion, the form of the Reynolds number to be 
used must be determined since the particle must traverse regions of variable unit 
Reynolds number. Since the test gas is first set in motion by the shock Sq and is then 
processed by an unsteady expansion, no uniform axial free -stream properties are avail- 
able with which to form a Reynolds number. The concept of an integrated Reynolds num- 
ber has been used in boundary -layer transition studies on cones (Nagamatsu and Sheer 
(ref. 8)) and on nozzle walls. A total integrated Reynolds number is defined as 


Rt= £ Rl 


where 

R i = J 

x P,i 


The total Reynolds number is based on the Reynolds number accumulated by a 
particle (as shown in fig. 11) as it traverses the region (2) (a region of constant free- 
stream properties), the unsteady expansion (a region of varying free -stream properties), 
and region (5) (a region of constant free -stream properties). Thus, Rj. takes the form 


~ R l, 



dx + 


' x p,(D 

x P ,f 


Rj(x) dx + 



Determination of the integrated Reynolds numbers in regions (2) and © is straight- 
forward. The procedure for carrying out the integration over the unsteady expansion fan 
was accomplished as follows. By use of the flow properties in region (2) as the initial 
conditions, thermodynamic properties based on frozen or equilibrium air were deter- 
mined through the expansion fan for 30.48 m/sec increments in the velocity until the value 
of U(§) was reached. Since the expansion fan is time -similar (conical) with respect to 
its origin, it is convenient to locate a new t',x f coordinate system, as shown in figure 11, 
with the origin at the apex of the centered expansion fan. Also, because the flow is conical, 
it is only necessary to determine one value of the integrated Reynolds number through 
the fan since the Reynolds number of any other particle traversing the fan can be scaled 
to that calculated value. Therefore, an arbitrary reference time at which a particle 
enters the fan tQ 6n t is picked and Xo >en t is determined from 

x o,ent = tojent^d) " a ©) 
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A Runge-Kutta type of integration is then carried out so that the position of the particle 
on each of the incremental U - a characteristics is known. This integration then 
gives the thermodynamic state of the particle as a function of its position and the final 
integration, 


R = 



R^(x) dx 


is done by Simpson's rule. Thus, the integrated Reynolds number of any particle travers- 
ing the fan can be characterized by this reference Reynolds number (Rf Q ), the reference 
times (t^ent and to, ex) at which the reference particle enters and exits the expansion 
fan, and the time (tg n ^ at which the particle of interest enters the fan. 

The integrated Reynolds number in region (D and expansion fan can be combined 
to give an analytic expression for the integrated Reynolds number at the trailing edge 
of the expansion fan. By reference to figure 11, the path of the shock Sq is 



( 1 ) 


and that of the particle in region @ is 


t’ 



+ b 


(2a) 


Also, 

v, _ x o,ent 

u o _ l o,ent " “Tf 

U © 

For some arbitrary value of te n t, 

v. *ent . 

b = 7i b o 

l o,ent 

and equation (2) becomes 

t' =YT- + ^- b o 
^ent 


(2b) 


Equations (1) and (2b) are solved simultaneously to give the point at which the particle 
was set in motion 
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o 


(3) 



u s,d) U (D 


The point at which the particle enters the fan is 

x p,ent = tented) - a (g)) (4) 

Subtracting equation (3) from equation (4) gives Ax pj the total distance that the par- 
ticle has traveled in region (2). 

The integrated Reynolds number in region (2) is simply 


R (2) = R f,(2) ( 5 ) 

Since the flow in the expansion fan is conical, the integrated Reynolds number through the 
fan scales is 


^ _ x ex ~ x ent 
Rf.o x o,ex " x o,ent 

where 

x ex = ( u © 

and 

x ent = ( u d) “ a (|))tent 


- a ©)f 


tent 


t; 


o,ent 


o,ex 


so that 


Rf = 


( u ® • a ®) zrB - ( u @ ■ a @) 


c o,ent 


tent 


R f,c t 

c o,ex " x o,ent 


( 6 ) 


(7) 


Therefore, the total integrated Reynolds number to the point at which the particle exits 
the expansion fan is 
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0 6X 

R ©,fan - < ( U ® - ' ( U @ " a ®| 


% 

1 1 
x o,ex _ x o,ent 


( u ©- a ©)-^fur 


-1 


o, ent \ U S, (I) U @ 


R l,©( l ent 


( 8 ) 


and that point is given by 

tpv = te nt 


L ex - rr 

l o,ent 


and 


x ex - ( u (g) " a (B)) t ex 


( 9 ) 


Since lines of constant Reynolds number will be of interest, they may be plotted 
by inverting equation (8) and solving for tg n t while holding Rt constant and marching 
through the expansion fan in increments of velocity so that Ut - at represents the 
termination of the expansion fan. Then 


tent.i - r r 


Rt 


A; 




Hf 


,0 


x o,ex - x o,ent 


( U @ - a ©) ' ~ 


to ’ ent ( u s!® _u © 


( 10 ) 


R 


■l,® 


and 






s I o,ex,i^, 

l ex,i ~ 7» r l ent,i 

\ to, ent/ 

x ex,i = (Ui ~ a i)tex,i 


(ID 


Constant Reynolds numbers in region 

/R - R/ 


will be on the line 


x’ = Xg x + 


(2),fan\ 


V 


(Rt > R @,tan) < 12 > 
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and 


t' 



(13) 


This procedure was carried out for the same conditions as used for figure 10. A 
series of constant Reynolds number curves were generated and the curve which seemed 
to best fit the data was Rf r = 40 x 10® as shown in figure 12. There is good agreement 
at the end points, but there is considerable difference between the data and the calculated 
Reynolds number about 9 meters downstream of the diaphragm. Haggard (ref. 5) has 
shown that for a range of operating conditions which includes this particular example, 
the early part of the flow processed by the expansion lags the local equilibrium state. 
Succeeding sections of the flow begin to approach, but do not reach, the equilibrium state. 
Consequently, in an attempt to estimate the effects of nonequilibrium flow on Rj., lines 
of constant integrated Reynolds number were generated by assuming that the flow was 
frozen at the region (2) condition. The Reynolds number 40 x 10 6 curve for the frozen 
case is plotted along with the previous data in figure 13. It may be assumed that the 
cases for relaxing flow lie between these two limits. Here, as with subsequent plots 
involving values calculated from the frozen state, only the equilibrium head and tail of 
the expansion fan are plotted, although there are small differences between the equilib- 
rium and frozen characteristics. The data points lie first on the curve determined by 
assuming frozen flow, and then shift to the curve representing the same total Reynolds 
number, but were calculated under the assumption that the flow is in equilibrium. 

To check the validity of the frozen to equilibrium transition and to check on unit 
Reynolds number effect, two other cases were run. The first used a value of pq that 
would produce a flow closely approximating the totally frozen flow. The results of this 
case are shown in figure 14. There is excellent agreement between a frozen transition 
Reynolds number of 10 x 10® and the experimental data. A second case, shown in fig- 
ure 15, used a value of pq which will approximate the flow in equilibrium. Again, the 
data points agree very well with a value of Rt r of 40 x 10®. However, since for this 
case there is little difference in the equilibrium and frozen states of the gas, it is dif- 
ficult to tell whether the flow is relaxing. 

For each run made, wall boundary -layer transition along the tube was found to 
occur at a single transition Reynolds number which was calculated by use of an integral 
technique. The integrated Reynolds number works equally well when the flow is in either 
the near -equilibrium or near -frozen flow regimes as well as when the flow transits 
between these two limiting cases. Thus, it is only necessary to determine one experi- 
mental transition data point in order to calculate the location of transition throughout 
the length of the expansion -tube acceleration chamber provided that point is determined 
before transition reaches the interface. 
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In all three cases shown, the constant Reynolds number curves representing tran- 
sition intersect the path of the interface and terminate at that point. This condition is 
due to a discontinuity in at the interface. There was no experimental evidence 
to indicate that the point of transition moved across the interface into the acceleration 
gas. In fact, the data show (fig. 15) that once the point of transition reached the inter- 
face, it continued to move with the interface. 

There are two possible explanations for this behavior of the transition point. First, 
had it been possible to make measurements much further downstream, the transition point 
might have moved across the interface. Because ^ » R^ it would take a parti- 
cle initially at rest and just downstream of the secondary diaphragm much longer to reach 
the transition Reynolds number than it would a particle just upstream of the diaphragm. 
Second, the structure of the wall boundary layers in regions @ and © is fundamentally 
different, one being made up of hot acceleration gas and the other consisting of an upper 
layer of cool test gas and a lower layer of hot acceleration gas. It is not realistic to 
expect one Reynolds number to determine transition in these two different regions, nor 
would one expect the point of transition to move continuously between boundary layers 
having such radically different structure. 

Since the path of transition does eventually intersect the interface, transition limits 
the available test time in an expansion tube, the maximum test time available occurring 
at the point where transition crosses the tail of the unsteady expansion. 

Some care must be taken when relating these data to the flow in similar facilities. 
The Langley pilot model expansion tube had unmachined, extruded tube walls which might 
affect the location of transition when compared with transition in a facility where care 
had been taken to maintain smooth tube walls. Also, the disruptive effects of transition 
on the free -stream flow, namely, the pitot-pressure dip, are related to the ratio of wall 
boundary -layer thickness to the tube radius. Consequently, if the conditions presented 
here were used in a tube having a much larger diameter than the tube used in this study, 
one would expect to see the available test time increase for the same running conditions 
because of the increased time required for the boundary layer to fill the tube. 

The effect of unit Reynolds number on the transition Reynolds number is shown 
in figure 16, where Rj r is plotted against R^ (gj. There is, as can be seen, some 
variation; however, the evidence is inconclusive because of the limited data available. 

It is of interest to note that although the spread in Rj. r is relatively small (lO x 10® to 
40 x 10®), these values bracket a very large area on the x,t diagram as illustrated in 
figure 14. 
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CONCLUDING REMARKS 


Wall boundary -layer transition has been measured experimentally in the Langley 
pilot model expansion tube. It has been shown that boundary -layer transition occurs 
along the length of the expansion-tube acceleration chamber at a constant transition 
Reynolds number up to the time at which the path of transition intersects the interface. 
The integral method of calculating the transition Reynolds number worked equally well 
for equilibrium flow, frozen flow, and flow in transition between these two limiting cases. 
Although transition criteria established in this paper may not be applicable to similar 
facilities, it has been established that if an experimental transition data point is given, 
the integral method developed in this paper will predict the location of transition along 
the tube. 

Transition of the wall boundary layer was found to be a factor in limiting the 
available test time in the expansion tube under certain operating conditions. However, 
it would be inappropriate to apply the present results directly to tubes of different diam- 
eter and wall roughness. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., February 5, 1974. 
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Figure 1.- Schematic diagram of expansion -tube flow cycle. 
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Figure 6.- Wall static pressure and radial distribution 
of pitot pressure. 
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= 8.274 kN/m 



Figure 7.- Center-line pitot pressure and pitot -pressure profiles at x = 6.035 meters downstream 

of secondary diaphragm. 
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of secondary diaphragm. 
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Figure 10.- Transition data as a function of particle path, 







Figure 12.- Transition data as a function of integrated equilibrium Reynolds number. 




Figure 13.- Transition data as a function of integrated equilibrium and frozen Reynolds number 





Figure 14.- Transition data as a function of integrated frozen and equilibrium Reynolds number 

Near -frozen case. 
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Figure 15.- Transition data as a function of integrated frozen and equilibrium Reynolds number 

Near -equilibrium case. 
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